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Abstract 



We calculate the Sivers distribution functions induced by the final-state interaction 
due to one-gluon exchange in diquark models of nucleon structure, treating the cases 
of scalar and axial-vector diquarks with both dipole and Gaussian form factors. We 
use these distribution functions to calculate the Sivers single-spin asymmetries for 
inclusive pion and kaon production in deep-inelastic scattering. We compare our 
calculations with the results of HERMES and COMPASS, finding good agreement 
for vr^ production at HERMES, and qualitative agreement for 7r° and production. 
Our predictions for pion and kaon production at COMPASS could be probed with 
increased statistics. The successful comparison of our calculations with the HERMES 
data constitutes prima facie evidence that the quarks in the nucleon have some orbital 
angular momentum in the infinite- momentum frame. 
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1 Introduction 



It is well-known that transverse-momentum-dependent distribution and fragmenta- 
tion functions, nowadays commonly referred to as TMDs, can have a nontrivial spin 
dependences and that the so-called "T-odd" TMDs can lead to single-spin asymme- 
tries [H |2[ 131 H] . They are also often referred to as "naively T-odd" , because the 
appearance of these functions does not imply a violation of time-reversal invariance, 
since they can arise through final-state interactions. The Sivers distribution function 
/fi^, schematically depicted in Fig. (H is the oldest example of such a function. 




Figure 1: Schematic depiction of the Sivers distribution function f^. The spin vector 
St of the nucleon points out of and into the page, respectively, and kx is the transverse 
momentum of the extracted quark. 

It describes the difference between the momentum distributions of quarks in- 
side protons transversely polarized in opposite directions. The Sivers effect was put 
forward as a possible explanation for the large single-spin asymmetries observed in 
p TT X experiments [HEIE]- Furthermore, it generates single-spin asymmetries 
in inclusive hadron production in deep-inelastic scattering (SIDIS) IHE]. These have 
been measured by the HERMES collaboration to be nonzero for Ip^ — /' vr X [7] , and 
updated HERMES results for pion and kaon production have been reported in [8]. 
More recently, COMPASS has analyzed their SIDIS data on pion and kaon production 
off a deutron target [S], and collected data using a proton target in 2007. Phenomeno- 
logical analyses of SIDIS data have been performed in Refs. [IHl [HI 1121 113] • The Sivers 
effect may also result, e.g., in asymmetric di-jet correlations 

inpTp ^ ietietX [TUflS], 
but these are not yet visible in the data analyzed to date [T6] . 

In recent years the importance in hadron physics of the role of the transverse 
momenta of the partons has been better recognized, since they provide time-odd 
distribution and fragmentation functions, and make possible single-spin asymmetries 
(SSA) in hadronic processes [U [6], [17]. Specifically, it has been understood that 
one-gluon-exchange final-state interactions (FSI) are a calculable mechanism for gen- 
erating a transverse single-spin asymmetry in SIDIS [6]. This FSI generates a Sivers 
effect when the distribution functions are allowed to be functions of the transverse 
momenta of the partons, as well as their longitudinal momentum fractions. Therefore, 
taking the transverse momenta of the partons into consideration enlarges the realm 
of investigation of the nucleon structure. 
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A simple scalar diquark model was used in [6] to demonstrate explicitly that this 
FSI can indeed give rise to a leading-twist transverse SSA in SIDIS, which emerged 
from interference between spin- dependent amplitudes for different nucleon spin states. 
It was observed in [6l [IE] that the same overlap integrals between light-cone wave- 
functions that describe the contribution to the nucleon anomalous magnetic moment 
from a given quark flavor also appear in the Sivers distribution for that quark flavor 
(with additional pieces in the integrand). Since these integrals are the overlaps be- 
tween light-cone wavefunctions whose orbital angular momenta differ by AL^ = ±1, 
non-zero orbital angular momenta of the quarks inside the proton are essential for 
the existence of the Sivers asymmetry [6l [18] . 

In this paper we calculate the Sivers distribution functions in SIDIS induced by 
the one-gluon exchange final-state interaction for models of the nucleon in which 
the spectator diquark is treated as either a scalar or an axial-vector. As we discuss 
below, the simplest SU(6) wavefunction of the nucleon suggests that the spectator 
diquark would be in a combination of these two states. We consider two types of 
form factor at the nucleon-quark-diquark vertices: the dipole form factor used in [19] 
and a Gaussian form factor. When either of these form factors is used, we find that 
at larger transverse momenta the asymmetry calculated using one-gluon exchange 
may exceed unity in magnitude, indicating that unknown higher-order effects must 
become important there □. Imposing the physical restriction that the asymmetry 
calculated by one-gluon exchange not exceed unity changes rather little the Sivers 
single-spin asymmetry after integration over the transverse momentum, indicating 
that the results we present should be quite reliable. The similarity of our predictions 
for dipole and Gaussian form factors also indicates the stability of our results. 

We compare our results for the Sivers single-spin asymmetries for 7r+, n^, vr^, 
and K- with the SIDIS measurements made by HERMES and COMPASS. We find 
good agreement with the HERMES results for 7r+ production, and qualitative agree- 
ment for 7r° and production. The experimental errors in the HERMES mea- 
surements of the Sivers asymmetries for vr" and K~ production do not permit any 
definite conclusions to be drawn. The same is true of the current measurements by 
COMPASS, and we look forward to increased statistics that could further test our 
predictions for pion and kaon production in SIDIS. However, the successful com- 
parison of our calculations with the HERMES data already constitutes prima facie 
evidence that the longitudinal projection of the net quark angular momentum in the 
infinite-momentum frame is non-zero. 



* A more complete description of the FSI can be made by introducing an appropriate Wilson- line 
phase factor in the definition of the distribution functions of quarks in the nucleon [TJ [501 [211 [HI [13] ■ 
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2 Basic Formulae for the Sivers Asymmetry 



2.1 Definition of the Sivers Asymmetry A^^^'^'^ '^^\x) 

The SIDIS cross section l + N^ ^ l' + h+ x on a transversely-polarized target contains 
8 spin-dependent azimuthal modulations. Here we consider only one of them - the 
so-called Sivers asymmetry. The relevant angular distribution of the cross section 
contains unpolarized (U) and Sivers (Siv) parts: 

(fa{x,y,z,PTAhAs) ^ d^(^uix,y,z,PT, (j)h) _^ d^asiv{x,y, z, Pt) ^-^/^^ _ ^\ 
dxdydzd'^PTd(ps dxdydzd'^PTd(f)s dxdydzd'^PTd(ps ' 

(1) 

where (ph and (ps are the azimuthal angles of the transverse momentum of the produced 
hadron and the transverse spin of the target relative to the virtual photon direction. 
The Sivers asymmetry is usually defined as 

.sin(<^^-0s)/ Id(j)hId(j)sd^a{x,y,z,PT,(l)h,(l)s)sin{(t)h-(j)s) 

St j d(f)hj d(f)sd^a{x,y,z,PT,(ph,(Ps) 

and we note that the integration singles out the au component in the denominator 
and a Siv in the numerator. 

Within the LO QCD parton model, we have 

(4) 



dxdydzd'^PTd(f)s 
where 



fi{x,'k±) and f^{x,'k±) are the unpolarized and Sivers quark distribution functions 
inside the nucleon, 0^ is the azimuthal angle of the active quark q, Dig{z,p±) is the 
unpolarized fragmentation function of q into the hadron h, and 

P± = Pt - zk± (6) 

is the transverse momentum of the produced hadron with respect to the fragmenting 
quark momentum. 

In order to obtain the dependence of Sivers asymmetry on any single kinematic 
variable such as x, one must integrate the unpolarized and polarization-dependent 
parts of the cross sections over the other three kinematic variables: 



Aa{x) = Int 



k, P 



(r(x) = Int 



M|Pt. 

C(a;, y)S,eJ | rf^k^ (x, k^) P^ - ^k^) 



where Int [...] denotes the following integration: 



Int 



G 



Zmax /--Pt max /"Stt 

rfz / ^IPtIIPtI / ' 

mill J Pt min J 



fcxi 



(i|k_L| |kj 



G 



(9) 



(10) 



The derivation of (IHl) is given in the Subsection 12.21 

In the later comparisons with experimental data, we use the following integration 
limits for asymmetries on a proton target at HERMES: 0.2 < z < 0.8, 0.05 GeV 
< |Pt| < 1.2 GeV, and for asymmetries on a deuteron target at COMPASS we use 
0.2 < z < 0.8, 0.1 GeV < |Pt| < 2.0 GeV, corresponding to the kinematic conditions 
of these experiments. 



2.2 Intrinsic Integration 

Let us consider two integrals: 

Ri = J /(kl)D(pl) 



R2 



■sm{<P^-<Ps) /(ki)D(pi). 



'Ill 



'12) 



The integrand of -Ri is a scalar function of vectors k^ and Pt in the two-dimensional 
transverse momentum space. This means that -Ri can be a function only of z and 
Pf,: 

i?i = ri(z,P|). 
On the other hand, f[T^ can be represented as 



i?2 — S1F2 — S2F1 



(13) 



Where are the components of the two-dimensional transverse spin vector S and 
the Fi are the components of 



/ d^k^ ^ /(ki)D(pi). 



(14) 



Since F is a two-dimensional vector, and the only vector remaining after integration 
in (fT^ is Pt, we have 

F = Pt$(P|), (15) 

and hence 



$(P|) 



Pt F 

p2 



(16) 



From Eqs. f[T3l - [T6|) we finally obtain 
where 

r,{z, Pi) = J d'k^ /(kl)I^(pl). (17) 

2.3 Calculations of the Sivers Function in Diquark Models 




Figure 2: Leading- order diagram contributing to fi in a diquark model. 




Figure 3: Leading- order diagram contributing to the Sivers function in a diquark 
model, via a one-gluon exchange final- state interaction (FSI). 

In this subsection we repeat the derivation of the Sivers function in the scalar 
diquark model with a constant Yukawa vertex given in [HI [21], and later generalize 
the result to the scalar and axial-vector diquark models with non-trivial form factors 
g{k'^) at the nucleon-quark-diquark vertex. We obtain, from fl56|) in the Appendix of 
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this paper or from (22) of [2l], the following expression for the distribution function 
in the Yukawa model: 

From Eqs. (28) and (30) of [24], we also have 

fclAi(x,k^) = -/^ (xM + m)(l-x) 

Here we set to unity the parameter a in Eq. (30) of [21]. Using f|T8|) and f|T9l) . we can 
write 

/iV(x,kx) 6162 M{xM + m) 



/i(x, k_L) 47r (xM + m)2 + ki 
where i?(x, k_|_) is given by 



klR{x,\iA.) , (20) 



,X X ^ ^ Zl _1_ ^^^^('^ . J ^ (21) 

We calculate the RHS of fl^ by first taking the direction of k^ as the x direction 
and using k and / to denote |k_L| and respectively, so that 

k± = A;(1,0), 

1_L = /(cos0, sin0), 

{l±_ ~ k±y = /cos0 - k, 

(I_L — k_L) = (Zcos0 — k, /sin0), 

(1± - kx)^ = + k^ - 2klcos(f) . (22) 

Then, / in (1211) becomes 

f 1 jiJA /cos0 — k 



I = — / ldld(j) 

TX J 



{P + B){P + k'^ - 2klcos^ 



°° , ., 1 /-I r^TT Icosd) - k 







After some calculations, we find that the last factor in fl23p gives 

fzl rl^ /cos0-fc N _ r 2 ^ when / < k . . 

U Jo 1- ^coscl>) ~ 1 when I > k ^^^> 



Using (^^, the expression fl23l) becomes 



2ldl 



k Jo (P + B) 

where we use 1/k = k'^/'k^, since we take the direction of k_L as the x direction. From 
(EID and ([25]), we have 

fl(..k,)x^ = ^/;'<i(li)|j^. (26) 

Then, using R{x,h±) in (!26|) . we can obtain the Sivers function from the formula 
(120|) . The result is identical to the result in [6l [2^ . 

Eq. (126!) is for the Yukawa model, presented in the Appendix, in which the 
nucleon-quark-diquark vertex is gik"^) = 1. We can generalize fl26p for a general form 
factor glk"^) = g{x, k^) by using the following formula: 



5" - ki io ^^^^^ WTb) ■ 



/^(.,k,)x^^ = -^ l/dill)^i^. (27) 



Using formula (1271) . we now calculate i?(x, k^) for the dipole and Gaussian form 
factors at the proton-quark-diquark vertex, respectively, obtaining the results shown 
in am and (HI. 



2.4 Generalized Diquark Models 

Whilst the simplest possibility for the spectator diquark is the scalar case, even in 
the absence of orbital angular momentum in the nucleon rest frame, it may exist 
in a spin-one, axial-vector state. Indeed, this possibihty is realized in the simplest 
non-relativistic SU(6) wavefunction for the proton: 
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In each of the terms in fl28l) and fl29l) . we have exhibited the spin states of the spectator 
diquark, showing exphcitly that it is in a combination of scalar and axial-vector 
configurations in different spin states. 

If the proton does have a naive SU(6) wavefunction in its rest frame, the dis- 
tribution functions of the u and d quarks inside the proton, and ff, are given 
by 

fi = ^/r + Ifi , ft = fi , (30) 

and those inside the deuteron are given by 

/deul ~ /deul = fl + fl ■ (31) 

The same relations also hold for the Sivers distribution functions, as we discuss below 
in more detail. 

In the following, we allow for non-trivial form factors at the nucleon-quark-di quark 
vertices, in both the scalar (s) and axial-vector (a) diquark cases: 

T, = g.ik') , = ^T'^Ts , (32) 

Two specific models for the form factors gsik"^) and Qaik"^) are discussed in the fol- 
lowing Sections. 



3 Dipole Form Factor 
3.1 Calculations 

We first consider a dipole model for the nucleon-quark-diquark vertex: 

*«(*^) = \/^ (F^A^yi = ^ (ki • '''' 

where R = s for the scalar and R = a for the axial- vector diquark, we used [191 Elj 

- k'ix, kl) = + -^Ml - xM\ (34) 

1 — X 1 — X 

and Br is given by 

Br = {1- x)K^ + xMl - x{l - x)M'^ . (35) 
We then have the following distribution functions: 

h {X, k^j - ^ (ki + ' ^ ' 
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yielding the following wlien integrated over k 

ci?^ N_ 7r{l~ xf[2{xM + mY + Br 

Tlie normalization factors Nr are fixed from the condition 



= ' ' ' ' — ^ . (37) 







^ dxf^{x) = 1 . (3J 



In the present cases of dipole form factors, we obtain from ([2 

1 1 r^'- ..2^ 1 



ki^ kl + B^^ bJ B^ikl + B^) ■ ^^^^ 

Thus, we obtain finally the Sivers distribution functions for the dipole form factors 
as follows: 

f±R(, W 1 - « "^"2 iV^,(l - xfMjxM + m) 
{x,k^)-an^ 16^35^(ki + 5^)3 ' 

In deriving (HOl) . we use for the gauge- field coupling to the axial- vector diquark in 
Fig. [3] the simple form ie2g"^{{P — I) + {P — k)Y, which is equivalent, for each 
polarization state, to the gauge- field coupling to a scalar diquark [J. Following [6], we 
fix ^ = —Cpois, where Cp = |- 

The values of are given by the overlaps of the proton wave functions of positive 
and negative helicities.We find: 

as = l and = (41) 

Details of the derivation are presented in the Appendix, where the relations between 
the SU(6) and light-cone wavefunctions are discussed, as well as the relations between 
light-cone and Bjorken-Drell spinors. 

We use as ~ 0.3 and choose the following values for the parameters of the model 
studied in this Section: 

m = 0.36 GeV, M, = 0.6 GeV, Ma = 0.8 GeV, A = 0.65 GeV. (42) 

We assume a Gaussian transverse-momentum dependence for the fragmentation func- 
tions: 

1 ^ 

^?,(^,pi) = — e^/^f,(^), (43) 



t We note that the resuhs for the axial- vector diquark in (j36|l . (f37|) and (|40|) are different from 
those of '^], which are obtained if J2\ = — is used. Ref. [26] considers various different 

possibihties for the polarization sum of the axial-vector diquark. More general forms of gauge-field 
coupling to the axial- vector diquark were used in Refs. [551 [23 127] ■ 
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with = 0.2 GeV^ as obtained in [28]. We use the leading-order fragmentation func- 
tions of [29] for the integrated fragmentation functions D^{z) . The integration 
is performed as described in Section 12.21 

We present in the top panel of Fig. H] the dependence of the calculated asymmetry 
on the upper limit in the integral. As one can see, the saturation of integral takes 
place around 2 GeV, which is a rather high value, though about 90 % of the integral 
is provided by the region \k±\ < 1 GeV. For definiteness, in this Section using dipole 
form factors, we use |kx|max = 2.5 GeV for the upper limit of the k^ integration. The 
leading-order perturbative QCD approach is believed to be applicable to SIDIS when 
all transverse momenta are much smaller than the virtuality of the hard scattering, 
Q. However, in SIDIS at fixed energy there is a strong correlation between the 
mean values of x and Q'^: for example, at HERMES {Q'^){x = 0.18) ^ 4 GeV and 
{Q'^){x = 0.28) ~ 6 GeV. Thus the highest virtuality at HERMES corresponding to 
the last populated x bins is of the same order of magnitude as the saturation value 
for the intrinsic transverse momentum integration. This consideration, together with 
the fact that the simple diquark model used here treats only the valence quarks, 
shows that this approach cannot be considered reliable for low values of x. Hence, 
our results should be considered as applicable only to x > 0.1. 

Another issue bearing on the accuracy of our results is that unitarity requires 
the analyzing power of the Sivers function to be less than unity in modulus for all 
values of x and k_L: Au^d = \k±f^''^{x, k±)/Mfi''^{x, k_L)| < 1. This is not always the 
case for the simplified leading-order one-gluon-exchange FSI that we consider, which 
tends to yield larger values at large x and k^, as seen in the middle panel of Fig. H] 
for scattering off a -u quark, and in the bottom panel for scattering off a c? quark □. 
The issue would be resolved if a full higher-order calculation were performed, but 
this is currently not available, so one must estimate the error incurred by including 
unphysical values in the integration over k^. In fact, we find that the naive predictions 
obtained by ignoring the unitarity issue differ from the modified calculations, that are 
obtained by truncating the k^ integration at the x-dependent unitarity limit shown 
in Fig. m by far less than the present experimental uncertainties, so this theoretical 
error may be neglected for the time being. 

3.2 Comparisons with Data 

In Figs. E] and M we compare the Sivers asymmetries for pion and kaon production 
calculated with a dipole form factor using = 0.3 and A = 0.65 GeV with data from 
the HERMES Collaboration obtained using a proton target. We see that our predic- 
tions for the TT"*" asymmetry (top panel of Fig. [5]) agree very well with the HERMES 
data, which exhibit a relatively significant positive asymmetry. The HERMES data 

related issue is that the asymmetry we calculate is large close to the contour in the (a;,kj^) 
beyond which unitarity is violated. 
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Figure 4: In the top panel, we show the Sivers asymmetry for tt'^ production off proton 
target for HERMES kinematics at x = 0.276 as a function of the upper limit of the 
|k_L| integration. The analyzing power of the Sivers function for scattering off a u 
quark (middle panel) and a d quark (bottom panel) as functions of x and = |kj_|. 
The calculated values exceed the unitarity limit in the white regions at larger x and 
kr- All these plots are obtained using dipole form factors. 
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Figure 5: Comparisons of our predictions for the Sivers asymmetries in the production 
of iT~^ (top panel), 7t~ (middle panel) and tt^ (bottom panel) with HERMES data, 
assuming = 0.3 and A = 0.65 GeV for the dipole model parameters. 



13 



-e- 


0.16 


!c 




-e- 


0.14 


.E i~ 








< 


0.12 




0.1 




0.08 




0.06 




0.04 




0.02 







In 

■o- 


0.2 










c 1- 

lA 3 


0.15 


< 






0.1 




0.05 









-0.05 




-0.1 






K" 


















1 


1 








1 


1 


1 


, 1 

1 1 1 1 


1 

III 


1 1 1 1 1 1 1 1 



0.05 



0.1 



0.15 



0.2 



0.25 



0.3 
X 



Figure 6: Comparisons of our predictions for the Sivers asymmetries in the production 
of (top panel) and K~ (bottom panel) with HERMES data, assuming = 0.3 
and A = 0.65 GeV for the dipole model parameters. 
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Figure 7: Comparisons of our predictions for the Sivers asymmetries in the produc- 
tion of 7T^ (top panel) and vr^ (bottom panel) with COMPASS deuteron target data, 
assuming as = 0.3 and A = 0.65 CeV for the dipole model parameters. 
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Figure 8: Comparisons of our predictions for the Sivers asymmetries in the produc- 
tion of (top panel), K~ (middle panel) and (bottom panel) with COMPASS 
deuteron target data, assuming ag = 0.3 and A = 0.65 GeV for the dipole model 
parameters. 
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for the 7r~ asymmetry (middle panel) are more equivocal, though they are compatible 
with positive values that are smaller than for tt~^, as predicted by our calculations. 
We also predict a positive asymmetry for 7r° production, which is in qualitative agree- 
ment with the HERMES data, as seen in the bottom panel of Fig. [5l In the case of 
the HERMES data for kaons, we see qualitative agreement for the asymmetry 
shown in the top panel of Fig. [6], though the experimental values are somewhat larger 
than the predictions, albeit with large errors. This difference could be explained by 
the fact that we do not include the s contribution, which dominates in fragmentation 
to the K~^. In the case of the K~ asymmetry shown in the lower panel of Fig. El the 
experimental values are similar to the predictions, though again with relatively large 
errors. 

The corresponding comparisons for COMPASS data on pion and kaon production 
are shown in Figs. [7| and [H respectively. In this case, again using a dipole form factor 
and as = 0.3 and A = 0.65 GeV, we predict small asymmetries for all three charge 
states for both pions and kaons. These predictions are compatible with the data 
available from COMPASS for 7r^,if^ and production. 

We conclude that our model gives qualitatively successful predictions for the Sivers 
asymmetries if a dipole form factor is assumed. 

4 Gaussian Form Factor 
4.1 Calculations 

In order to assist in evaluating the reliability of our results, in this Section we make 
a similar analysis of the predictions for the Sivers asymmetries obtained using a 
Gaussian form factor at the proton-quark-diquark vertex given by 



where —k'^ = -H^ + i^Mn — xM"^ . In this case, we have the distribution functions 




(44) 




/o(x) [{xM + mf + kl 




A?{1-^) 




(45) 



Al{l-x) 



where 



x(l-x)M^ -xaA 



fo{x) = e 



(46) 



After integration over k_L, we obtain the distribution functions 



(x) = TiNn /o(x) [(xM + mf + Aj 



(1-^)] , 



(47) 
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where the normahzation factors Nr are fixed by 

Cdxf^{x) = l. (48) 
Jo 

We obtain from (1271) for the present Gaussian form factor 

R{x, kx) X e^^^^ = TT I ^(li) e"^^^?(^ = 2A?(1 - a;) (l - e'^^^A . 

k]_ JO kj_ ^ ^ 

(49) 

Finally, we obtain the following Sivers distribution functions for the Gaussian form 
factor: 

/iL/(a:,k^) =aR^NR 2M{xM + m) fo{x) e~^^ ( 1 - e~^^ ) , 

(50) 

where = 1 and = — | as in fHOl) . Relevant formulae and details of the derivation 
and the relation between the light-cone and rest-frame SU(6) descriptions of the 
nucleon wavefunction are given in the Appendix. 

As one can see in the top panel of Fig. [9], the convergence of the k^ integration 
is faster with the Gaussian form factor than with the dipole form factor, and in this 
section we use |kx|max = 1-5 GeV for the upper limit of the k^ integration. In 
the middle panel of Fig. [9l we delineate the regions of the (x, k^) where the naively 
calculated one-gluon-exchange FSI contribution to the Sivers asymmetry in scattering 
off a M quark exceeds unity (white), and also indicate the magnitudes calculated at 
lower values of x and k^. The bottom panel presents the same information for 
a d quark. As in the case of the dipole form factor the unitarity cut of the k^ 
integration limit alters our predictions by amounts that are again far smaller than 
the experimental errors. 

The bottom panel presents the same information for a d quark. As in the case of 
the dipole form factor, we cut off the k^ integral at the solid line, which alters our 
predictions by amounts that are again far smaller than the experimental errors. 

4.2 Comparisons with Data 

We display in Figs. [10] and [11] comparisons of our predictions for the Sivers asymme- 
tries with HERMES data on pion and kaon production, respectively. These predic- 
tions are calculated with the Gaussian form factor using as = 0.3, Ai = 0.5 GeV. As 
in the dipole case, we find quantitative success for the 7r+ asymmetry and qualitative 
success for the tt" asymmetry, while the experimental errors in the vr" asymmetry do 
not permit a firm conclusion to be drawn. In the case of the asymmetries, we 
again find a qualitative success for the case, though the measured asymmetry 
is larger than our prediction. This may again be explained by ignorance of the s 
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Figure 9: In the top panel, we show the Sivers asymmetry for vr"*" production off a 
proton target for HERMES kinematics at x = 0.276 as a function of the upper limit 
of the |k_|_| integration. The analyzing power of Sivers function for scattering off' a u 
quark (middle panel) and a d quark (bottom panel)as functions of x and kr = \^±\- 
The values exceed the unitarity limit in the white regions at larger x and kr- All these 
plots are obtained using Gaussian form factors with parameters Ai = 0.5 GeV and 
as = 0.3. 
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Figure 10: Comparisons of our predictions for the Sivers asymmetries in the produc- 
tion of (top panel), n~ (middle panel) and (bottom panel) with HERMES data, 
assuming as = 0.3 and Ai = 0.5 GeV for the Gaussian model parameters. 
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Figure 11: Comparisons of our predictions for the Sivers asymmetries in the pro- 
duction of (top panel) and K~ (bottom panel) with HERMES data, assuming 
as = 0.3 and Ai = 0.5 GeV for the Gaussian model parameters. 
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Figure 12: Comparisons of our predictions for the Sivers asymmetries in the pro- 
duction of (top panel), tt~ (bottom panel) with COMPASS deuteron data for vr^, 
assuming = 0.3 and Ai = 0.5 CeV for the Caussian model parameters. 
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Figure 13: Comparisons of our predictions for the Sivers asymmetries in the produc- 
tion of (top panel), K~ (middle panel) and (bottom panel) with COMPASS 
deuteron data, assuming = 0.3 and Ai = 0.5 GeV for the Gaussian model param- 
eters. 
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contribution. In the K~ case, the prediction is also quahtatively successful, though 
no definite conclusion can be drawn. 

The corresponding comparisons between our predictions and COMPASS data on 
pion and kaon production are shown in Figs. [12] and [T31 respectively. In the cases 
we predict smaller Sivers asymmetries than for HERMES, and the data certainly 
reflect this trend, though the data are quite compatible with zero. In the cases of the 
kaon asymmetries shown in Fig. [131 the predicted asymmetries are very close to zero, 
as are the values measured by COMPASS. 

We see that our model gives qualitatively successful predictions for the Sivers 
asymmetries also if a Gaussian form factor is assumed, and that the predictions 
using this and a dipole form factor are qualitatively similar. This gives some further 
confidence in the stability of our results and the conclusions we draw. 

5 Conclusions 

We have studied the Sivers single-spin asymmetry in SIDIS generated by the mech- 
anism of a one-gluon-exchange final-state interaction. We derived a general formula 
that can be used to calculate the Sivers distribution function for diquark models hav- 
ing different form factors at the nucleon-quark-diquark vertex. We calculated the 
Sivers distribution functions in diquark models with both dipole and Gaussian form 
factors, and compared the corresponding predictions for Sivers single-spin asymme- 
tries in pion and kaon production in SIDIS with the results of HERMES and COM- 
PASS. The predictions made using dipole and Gaussian form factors are quite similar, 
and are relatively insensitive to the unphysical values of the model calculations at large 
X and k^. 

We find qualitatively successful results for the asymmetries in h^^tt^ and 
production. In the case of production at HERMES, the measured values are 
even larger than our predictions, refiecting the possible importance of an s contribu- 
tion. In other cases, particularly at COMPASS energies, many of the experimental 
measurements are currently compatible with zero, and greater accuracy will be nec- 
essary to confront our theoretical predictions. On the theoretical side, it is desirable 
to improve the accuracy of our predictions, in particular by going beyond the sim- 
ple one-gluon-exchange final-state interaction. However, the relative success of this 
first confrontation between HERMES and COMPASS data and our naive predictions 
is an encouraging indication that one may be able to understand satisfactorily the 
Sivers asymmetries in SIDIS, which are rather subtle aspects of hadron dynamics in 
deep-inelastic scattering. 

One prima facie conclusion from the successful comparison between our predic- 
tions and the HERMES data on tt"*", 7r° and production is that the quark partons in 
the nucleon must have nonzero orbital angular momentum in the infinite-momentum 
frame. 
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Appendix: Wavefunctions of Scalar and Axial- vector 
Diquark Models 

The expansion of the proton state in terms of light-cone Fock states is 



X ipn{xi, k±i, Xi) n; XiP^ , x^P^ + kj_i, A, 

where the light-cone momentum fractions Xi = kf / P+ and k_Lj represent the momenta 
of the QCD constituents. The physical transverse momenta are p^j = XjPj, + k^j. 
The Xi label the light-cone spin projections of the quarks and gluons along the 
quantization direction z. The n-particle states are normalized as 

n 
i=l 

Here and in the following we do not display the other quantum numbers of the partons, 
i.e., color and quark flavor. 

In order to construct diquark models, we take the following form factors at the 
proton-quark-diquark vertex, for the scalar (s) and axial-vector (a) diquarks respec- 
tively: 

T, = g,{e) , = ^7^5 . (53) 

We can then obtain the light-cone wavefunctions of scalar and axial-vector diquark 
models from Fig. [TH In this Appendix we consider elementary vertices given by 
Qsik'^) = Qaik"^) = 1- In Sections [3] and m we introduce dipole and Gaussian form 
factors, respectively. 
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Figure 14: Diagram giving the light-cone wavefunctions of scalar and axial-vector 
diquark models. 



A. Scalar Diquark Model 

We also use the term 'Yukawa model' for the scalar diquark model described in this 
Subsection. Each Fock-state wavefunction of the physical proton with total spin 
projection = ±| is represented by a function ip^'^ {xi, k±i, Aj), where 



— {kf , k_|_j, fcj ) 



(54) 



specifies the momentum of each constituent and Aj specifies its light-cone helicity in 
the z direction. 

From Fig. [THwith the scalar vertex T^, the = +| two-particle Fock state in 
the scalar diquark model is given by 



v|/T ( p+ p 

two particle V ) - 
dx (i^k I 



x(l — x) IGtt^ 







ij]i{x, kj 



(55) 



where 



# i(x,kj 
ip^_i{x, kj 



(xM+m) 



(+fci+ifc^) 



(56) 



The scalar part of the wavefunction cp is given by [301 EI] 



(p{x, kj 



9 



1 



X 



X 



X 1— x 



5 



(57) 



where 



B = x 



(l-x)(-M^ + ^ + -^). 
^ X 1 — 



(5J 
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Similarly, the = —■^ two-particle Fock state is given by (SOj [3T 



i 

two particle 



(P+,Px = Oj 



(59) 



dx (i^k I 



where 



^ii(a;,kj 
^ii(a;,kj 



(-fcl+ifc2) 



(xM+m) 



--; xP+, kj 



(60) 



The quark-state bases on the right-hand sides of (l55l) and (l59l) correspond to light- 
cone spinors. The Bjorken-Drell spinors u^^{k) and the light-cone spinors u^^(/c) are 
related by 



<?(fc) 

^ 2 



1 



{k+ + m)2 + 
1 



(/t+ + m) M^<i(A:) - (A;^ + zA;2) u^:'i(A;) 



( -(-k^ + 1^) ufrik) + (k+ + m) ufUk)) . (61) 

V "'"2 2 / 



'(fc+ + m)2 + fc2 
The Bjorken-Drell spinors u^^{k) satisfy 

j^ulY{k) = 0, ru^Y{k) = 0, rulY{k)=u'tm, J^u^Yik) = ulYik) , (62) 

^2 2 ^2 2 2^2 

where = ± ij"^ and 



= + 



a' 

(T^ 



d 

dk^ 



6^23 = 1 . 



(63) 



Since k~^ = xM in the proton rest frame, from ( 1611) the Fock states given in (1551) and 
(13^ with the light-cone wavef unctions and (1^ correspond to +| ; a;P+ , k 



and 



2 i , kj 



BD 



BD 



, respectively, in the proton rest frame. 



B. Axial- vector Diquark Model 

For the polarization vectors of the axial-vector diquark appearing at the vertex of 
Fig. [TU we use the following set of three polarization vectors: 



(e(+l)o,e(+l)^e(+l)^e(+^)3) 



(0,-1,-2,0) 



(0,+l,-z,0) 



(64) 



pO 
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These polarization vectors are those used in Ref. [19] which satisfy Y.x = 
~g^iv + PfiPu/M"^- In the proton rest frame, e^^^'^ is the spin vector oriented in the 
z-direction and e'^^^^'^ are those circularly polarized in the x-y plane. 

From Fig. [14] with the scalar vertex ■ e*-^^*, the two-particle Fock state for the 
proton with = +| has six possible spin combinations for the quark and axial- vector 
diquark: 



two particle V 



0, 



+ V'|io(a;,kx) 
+ i^li_^ix, k± 



+- + 1 ; xP^ 
+ ^ -1; xP' 



x{l — x) IGir^ 
. kx^ + V^li+i(a;,k_L) 



(65) 



-- +1; xP' 



.) + i^[ioix, k^) 



— 0; xP^ 



1 



- - 1 ; xP+ 



where the two-particle states \sf,sl; x,'k±) are normalized as in (1521) . Here sf and 
denote the z-component of the spins of the constituent fermion and boson, re- 
spectively, and the variables x and k^ refer to the momentum of the fermion. The 
wavefunctions are given by 



2 

t ^ 

ip.iAx, k±] 



2 i-k^+ik^) 

3 X ^ 1 

2 (xM+m) 

1 (xM+m) 

3 X ^ ' 
r (+fci+»fc^) 
3 X 



(66) 



0, 
0, 



where 



V9(x, kj 



l-x 



(67) 



Similarly, the wavefunctions for a proton with negative helicity are given by 



2 

2 ^ 
^2 

'4)liAx,ki_] 



= + 
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1 i-k'+ik^) , 

3_^r^<^ ' 

2 (xM+m) 

3 X ' 
2^ (+A:^+ifc^) 
3 X 



(68) 



The coefficients of in (l66l) and (!68|) are tlie matrix elements of ^« ' 

^(A)* HP^^_^^^P±) Tf^i^io}^ aj-g ii^Q numerators of the wavefunctions corresponding to each 
constituent-spin configuration. 

Because of the transformation properties (!62|) of the Bjorken-Drell spinors u^^{k), 
when we construct the proton spin states \P; A = ±|) by using the Clebsch-Gordan 
coefficients with u^^{k) for spin-half components, they satisfy 

J+|P;A = +i)=0, J-\P;\ = -^)=0 

J-\P;X = +l) = \P;X = ~l) , J+|P;A=-^) = |P;A = +i) , (69) 

where J* is the angular momentum operator for the proton given by J* = J^aJa — 
J^ai^a + ^a)) wlfich is the sum over the constituents a, and J"^ = ± 

Since k'^ = xM in the proton rest frame, from (!6T!) the Fock states with the 
light-cone wavefunctions given in (!66l) and (!68l) correspond to the following states in 
the proton rest frame, respectively: 




The above states are angular momentum eigenstates |P;A = +|) and |P; A = — |), 
respectively, which satisfy (|^ . In the proton rest frame, the states described by 
(|66|) and ( fl68l) are identical to the axial- vector diquark states appearing in the SU(6) 
wavefunctions fl28l) and fl29|) . Therefore, in the proton rest frame the wavefunctions of 
the axial- vector diquark states derived here through the vertex given in Fig. [T^ with 
the polarization vectors given in flMl) coincide with those in the SU(6) wavefunction 
for the proton. 
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